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Abstract 

We calculate the large Yukawa coupling corrections to the top and bottom scalar 
quark pair production in e"^e~ annihilation within the Minimal Supersymmetric 
Standard Model. We include the vertex corrections and the corrections to the 
gauge boson propagator enhanced by large masses. We find the total corrections 
are quite significant. In some regions of the parameter space the corrections are 
larger than 10%. 

PACS numbers: ll.30.Pb, 14.80.Cp, 13.85. Qk, 12.15.Lk 



I. INTRODUCTION 



Supersymmetry(SUSY) is one of the most attractive extensions of the Standard Model(SM). 
It provides an elegant way to stabilize the huge hierarchy between the electroweak and the GUT 
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scales against radiative corrections |T]]. Moreover, supersymmetric models offer a natural solution 
to the Dark Matter problem |^ and allow for a consistent unification of the all known gauge 
coupling constants in contrast to the SM f^. Due to the theoretical appealing of SUSY, the 
search for supersymmetric particles is one of the main issues in the experimental programs at the 
CERN e^e^ collider LEP2 and Fermilab Tevatron It will play an even more important role 
at the future Large Hadron Collider (LHC) [^] and the Next e+e~ Linear Collider 

Although the colored supersymmetric particles, squarks and gluinos, can be searched for most 
efficiently at hadron colliders, for a precise determination of the underlying SUSY parameters 
lepton colliders will be necessary. For the experimental search it is useful to predict the production 
rates of these particles precisely incorporating radiative corrections. Up to now, many works have 
been devoted to the QCD corrections to various sparticle production rates. QCD corrections 
to colored sparticle (except stop) production at hadron colliders were discussed in detail by W. 
Beenakker et al. 0. The corresponding corrections to the top squark production were given in 
another paper ||^. The QCD and SUSY-QCD corrections to non-colored sparticle production at 
hadron colliders were given in 0] and those to colored sparticle production at e^e~ colliders were 



given in |[T0|. 



In this paper, we consider the electroweak corrections to the third generation diagonal squark 
pair production in e+e~ annihilation, e+e^ titi,bibi, due to large Yukawa couplings. Our 
framework is the Minimal Supersymmetric Standard Model(MSSM) As is well known that 
there are five physical Higgs bosons in the MSSM, two CP-even neutral Higgs bosons, one CP-odd 
neutral Higgs boson and a pair of charged Higgs bosons. Their supersymmetric partners, higgsinos, 
are components of two charginos and four neutralinos in the MSSM. The top and bottom squarks 
have Yukawa couplings with these Higgs bosons and higgsinos, which are proportional to rrit cot j3 
or rrihta.njS, where tan/3 = V2/V1 and Vi, V2 are the vacuum expectation values of the two Higgs 
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doublets. These interaction terms are large in the region of small or large tan/? and they can 
even be leading electroweak corrections for tan/5 ~ 1 or tan/? ~ mt/mf,. On the other hand, the 
internal gauge bosons may also have large corrections enhanced by large masses due to virtual 
heavy particle loops such as the top or stop loops. For consistency, we also include such corrections 
in our calculations. We calculate in the 't Hooft-Feynman gauge. We find the total corrections are 
quite large in some regions of the MSSM parameter space allowed by present experiments, which 



can be larger than the SUSY-QCD corrections to the same process due to gluino exchanges [10 



This paper is organized as follows. In Sec. II we present the renormalization scheme adopted in 
our calculation. Some analytic results are given in Sec. Ill and the numerical results are discussed 
in Sec IV. Finally we summarize the conclusion in Sec V. The relevant pieces of the Lagrangian 
are presented in Appendix A and some analytic expressions are collected in Appendix B. 



II. RENORMALIZATION SCHEME 

In this section we briefly discuss the renormalization scheme adopted in our calculations. To 
calculate the electroweak corrections to the process e~^e~ titi{bibi) at one loop level, we do not 
need to consider the renormalization of the Higgs sector after imposing the vanishing of the tadpole 
terms. (However, we adopt an approximate Higgs mass formula including radiative corrections.) 
Thus, the renormalization scheme is focused on the gauge sector. It differs only slightly from that 



given by M. Bohm Another complexity arises from the renormalization of the squark mixing 
angle. 
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A. Gauge boson renormalization 



The diagonal production of squark pairs proceeds through S-channal photon and Z boson ex- 
change at tree level(see Fig. 1). The longitudinal part of Z boson does not give any contribution to 
the process. In the MSSM the photon and Z boson may mix with the CP-odd neutral Higgs boson 



and the neutral Goldstone boson G"^ at one loop level [T^. However, for diagonal production 
of squark pairs, such mixing does not give any contribution either. So only the renormalization of 
the transverse part of the gauge bosons is needed. 

To respect gauge symmetry explicitly, each gauge multiplet is associated with one renormal- 
ization constant 

Wl - (Zf , B, {Zffl^ B, , 

g,^ZY{Z^r'l-g,^ g,^Zf{Zf)-^l-g, . (2.1) 



The Weinberg angle is defined by the on-shell condition cos^iy = ^^j, where and Mz are 
the masses of W and Z bosons. Now we denote 

ciy = cos , svK = sin6'vy (2.2) 



as abbreviations throughout the paper and 

= sl^bZf + ^y,bZf , 5Zf = s^^bZf + cl^bZf , 
(5Z7^ = -c^s^^bZ'f - 6Zf) (2.3) 

as the renormalization constants for the photon, Z boson and 7 — Z mixing terms respectively. 
Then we get 
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.A. 



1 + \6Zl -6Zf + 5Zf ^ 
5Zf - 25Zf 1 + UZl 



.A, 



from which we can see the 7 — Z mixing term. 



(2.4) 



The renormahzation constants Z^2 fixed by the following on-shell conditions 



S5r(M2,) = S|(M|) = SJ^ (0) = 0, 



(2.5) 
(2.6) 
(2.7) 



where the Sys represent the renormalized self-energies and the T"'^'^ represents the renormalized 



photon-electron vertex. From Mw = 5'2/2-\/f 1 + and Mz = hy gf + giy + V2 we get 



Ml 



■'W 



{26 Z^ - 36 Z^) - {26 Z^ - 36 Z. 



(2.8) 



Throughout this paper we shall keep only corrections proportional to a large mass M > Mz- 
All terms independent of the large masses or depending on them only logarithmically will be 
ommitted. It is found that no terms proportional to large masses enter the expressions ^'^^2 ^ \k^=o 
and J]'^^{0). The same is true for 6Z^ determined from (2.6). Taking into account these facts we 
obtain from the renormahzation conditions (2.5)-(2.7) 

S^(A;2) 



Sw6Z^ + (^6Z^ — 6Z2 — - 
S^^(O) 



-6Zf + 6Zl^ 

6Zl = 



A;2 



|fc2=0 











(2.9) 



The calculations here are similar to those in |T2|. The four equations in (|2.8| ) and ( p.9|) completely 
determine all wave function renormahzation constants as only four of them are independent. We 
get from these equations 



.yz ^ cw - sir ( 5M\ _ 5Ml, \ 



The self-energies of gauge bosons S^, S^^ and Sf (Fig. 2a) for A;^ ^ may contain terms 
proportional to large masses. However, it turns out that their contribution to renormalized gauge 
boson propagators can be neglected. As an example, let us look into the top-quark loop correction 
to the renormalized Z boson propagator. This propagator can be written as 

The contribution of the top-quark loop to Sf (A;^) is written down in (B.20). Although contains 
terms proportional to m^, it can be checked that the combination (^t(^^)^^t^(^'^z)) jg enhanced 



by and can be negelected compared to 6Z2 in ( p.ll|) for all values of k'^ (See discussion following 



(B.21)). Thus, Z and 7 — Z boson propagator can be written as 

:i-5Zf) (2.12) 



and 



^ bZf (2.13) 



respectively. 

The analytic expressions for the gauge boson mass corrections calculated from Fig. 3 are given 



in Appendix B. By using (B.15)-(B.19) we have checked that the divergences in individual terms 
of the expression (B.4) for ^ti¥- — -rp^ are cancelled out after omitting terms not enhanced by a 

z 

large mass M > Mz. Hence 5Z^ and 5Zf obtained from (p30| ) are finite. This is also confirmed 



numerically. It should be noted that the full expression for and SZ2^ are of course divergent. 
The finite results obtained here are consequences of omitting divergent terms independent of large 
masses M > Mz- After cancellation of divergences in the full expression such terms can not 
induce finite corrections proportional to a large mass. 



B. Renormalization of squark wave function 



There are two scalar partners ql and qr for every quark q in SUSY theories. They mix and 
form two mass eigenstates qi, q2 which are related to the original fields by 



ML 



(2.14) 



where 



(a ■ 

cos Oq — Sm Oq 



y sin Oq COS 9q I 



(2.15) 



We will adopt a scheme in which both stops and sbottoms are defined on shell. We give the 
formulas for stops here while those for sbottoms are similar. 

The complexity of the squark wave function renormalization is due to the fact that the two 
diagonalized states q\ and qi mix again at one loop level(See Fig. 4). Write the bare stop fields as 
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(2.16) 



(We use 5Zf and SZfj to represent the wave function renormalization constants and Zg the mixing 
matrix.) The on-shell renormalization conditions require that the mass parameters are the physical 
masses, the residues of the squark propagators on shell are one and the mixing between on-shell 
squarks should be absent, i.e. 



E"K_)=0 , E"(mi) = , 8=1,2 , 



d 



0. 



^ t2 



(2.17) 



From the above equations, we get 



Sml = E(m|) , 



6Z^ 



-S'(m|) , 



6Zi 



— 



(2.18) 



where S'(p ) is the derivative of S(p ) with respect to jr. 

The wave function renormahzation constant matrix can be decomposed into a symmetric and 
an antisymmetric part 



V2 



\{6Z\^ + 6Z\^) l + Uzl 



\ 



I 



I 



1 



\[bZ\^-bZ\^'^ 



V"2 



\[bZ\2 - 6Z21] 



1 



[2.19) 



where the off-diagonal elements of the symmetric part are ultraviolet finite and the antisymmetric 
part can be interpreted as a rotation matrix in the first order. Besides the wave function and 
gauge coupling constant renormahzation defined above, an additional renormahzation of the stop 
mixing angle 9i ^ 9i + 59i must be introduced to make the Ztitj vertex part finite beyond the 
tree level. We choose 59i such that this additional rotation just cancels the last factor in ( |2.19|) , 
that is. 



69i 



5Z\2-5Zi, 



(2.20) 



This is the same scheme as used in It is found that with this choice of mixing angle renor- 
mahzation the ultraviolet divergence in the vertex graph is exactly cancelled. 



The analytic expressions for the self energies S^-'s calculated from Fig. 4 are given in Appendix 



B. 



C. Renormalization of the gauge boson and squark vertex 



With the choice of (|2.2CI| ) the complete one-loop electroweak corrected Lagrangian for the gauge 
boson and stop interaction vertex is given by 



L --- 

Ztiti 



2 . <^e 1^^. 

,e 1 + - + 6Zr + -5Zl 



le 



SwCw 



(2.21) 



le 



SwCw 

ie 

SwCw 

2 



1 2 



e cos 6w sm 611^ 2 / 



- sing,- cos ^ /^^^t^^^^ 1 t,: UZ 



-ie\ 



-5Zf + 5Zf)i: UZ^ . 



(2.22) 



By using (|2.9|), the above two equations are reduced to 



C --- 

Ztiti 



■-ie{l + 6Z^Jii* 9^ M 



26 



1 9 



'l + 6Z^_)t7 d^'UZ^ 



sin 6*, cos 9^ 



5Z{^ + 



(2.23) 



(2.24) 



The corresponding Lagrangian for the sbottoms is similar. 



III. ANALYTIC RESULTS 
A. Vertex corrections 

All the quantum effects, including the vertex corrections and the corrections to the Z gauge 
boson propagator but not the 7 — Z mixing contributions, can be written in concise forms by 
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defining two effective coupling constants D'^ and . Let 



p7« 



■^Ztt 
eff 



-%e 



+ + 



le 



SwCw 

ie 

SwCw 



2 



(3.1) 



sin6'iCos6'i^ — i^— ^ + 



(3.2) 



wliere A''' and A^ are tlie vertex corrections by exclianging virtual Higgs bosons, cliarginos and 
neutralinos to the 7 and Z vertices respectively, as depicted in Fig. 5. It should be noted that 
the contribution from Fig. 5(f) is zero. The contributions coming from Fig. 5(d) and 5(e) cancel 
each other because the g — g — A" coupling changes signs when the momentum of the squark 
changes signs. Our results are analytically and numerically confirmed by that the UV divergence 
are cancelled precisely as it should be. The results have been verified by testing the Ward identity. 
The analytic expressions for A''' and A^ are listed in Appendix B. 



B. cross section 



We need to consider two types of contributions, one is the vertex corrections and the other is 
the corrections to the internal propagators shown in ( |2.12| ) and ( |2.13| ). We denote the amplitude 
due to the mixing as T'^^ and T^"^ where T'^^ corresponds to the photon propagator on the 
left and the Z boson propagator on the right in Fig. 2 and T'^^ the vice versa. We use T'^ and 

to represent the amplitudes calculated from ( p.lj) and (|3^ ). Then the cross section can be 
expressed as 



Tca 



Am 



2x1 



a 



L + M 



Mi 



(3.3) 



where s = {pi + P2) is the s-channal Mandelstam variable(Fig. 1) and L, M and are 
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n2 



5 



M = (1 - 4^2^ + 8^^ 



+ 



z\2 



'1 -4s 



2 



N = ^^(1 - 4s^) + , ^ M l - 44) + • 

S^C^ ZCw-Sw SwCw 

L comes from the square of T"^. The first term in the square bracket of M comes from the square 
of while the second term in it comes from the interference between T^"^ and . The last 
term in M comes from the interference between T"^^ and T^. The first term in is due to the 
interference between and T^, the second term is due to the interfernce between T^'^ and T'^ 
and the last term is due to the interference between T"^^ and T'^. 

Throwing away the one-loop corrections, we regain the tree-level formula. 



C. Higgs boson mass formula 



The Higgs sector is strongly constrained by supersymmetry |T5|. In the tree level a light Higgs 
boson exists with an upper mass bound Mz- Radiative corrections can considerably change the 
Higgs mass spectrum. In our calculations we adopt an approximate Higgs mass formula which 
incorporates the one loop radiative corrections. It is given by 

^^-^HO,hO,eff — ^ 



(M^ + M|)^ - M1,.M| cos= 20 + ^(Ml, - M|) (3.4) 



where 



N^Gprnt ( mi^nij^ A,{A, + /i cot /?) m\ 

^* ~ rR ■ 2 n 2 ^ 2 2 

V 27r2 sm'' (5 \ mi - m| m^-^ 
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^ A^At + fi cot / ml + m| m-,^ 
(m| - m|)2 \ ml - m| m^-^ 

Let cut = we return to the tree level formula of neutral Higgs boson masses. Although the 
correction about Higgs masses is a two-loop effect to the squark pair production, we find it can 
greatly affect the numerical results. 

The charged Higgs boson mass is given by 

M|± = Mlo + . (3.6) 




IV. NUMERICAL RESULTS 

Now we turn to discuss the numerical results. Since the cross section of squark pair production 
is sensitive to the squark masses, we use the two stop masses, m^-^ and m^-^, as input parameters. 
Making the following assumptions for simplicity 

A = A , (4.1) 

where the ms and As are the scalar masses and trilinear soft breaking parameters, we are left 
with only two free parameters in the squark sector. The sbottom masses are then determined by 
mi^, mi^, fi and tan/3. For simphcity we also assume the GUT relation mi = (5/3) tan^ 6'vi/^2 
where mi and m2 are U(l) and SU(2) gaugino masses respectively. The chargino and neutralino 
sectors are determined by taking m2 as another free parameter, m^o and tan/3 determine the 
MSSM Higgs sector. These free parameters are constrained by the experimental mass bounds. 
We impose m^^o > 90GeV 0, m^o > 35GeV, > 95GeV |18|,|T9[ and m^^ > 150GeV. To 
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discuss the large Yukawa couplings we focus our attention on the regions of small and large tan (5. 
The MSSM may seem unnatural for these values of tan/? |^^. However, they are actually not 
excluded by present experiments even for m/^o > 90G'eV^(See e.g. Ref. Other parameters 



are taken as a = 1/128, Mw = SOAGeV, Mz = 91.2GeV, mt = 174GeV, = 4:.7GeV and 
sin^ ^VK = 0.223. 

In Fig. 6, we show the cross section a{e^e^ titi,bibi) as a function of the collision energy 
^/s for = ISOGeV^, mj^ = 450GeV^ and fx = m^o = m2 = 4006*61^ for small and large tan/3 
scenarios. For tan/3 = 1.5 the two sbottom quarks are almost degenerate. However, for tan j3 = 30 
the lighter sbottom can be as light as ti and its production rates are much larger than those of 
the heavier one. 

We then calculate the corrections to the cross section a{e~^e~ titi) at ^/s = 206GeV at which 



LEP2 can run in 2000 |£2|. In Fig. 7, we show 5a /a as a function of the parameter /i by taking 



= 92GeV, which is slightly heavier than the present lower limit ||19|J23|1 , and m^^ = 3506*61^ 
for tan/3 = 1.5, my^o = 400, SOOGeV^ and 1712 = 200, SOOGeV^. We can see that the corrections are 
not sensitive to m2. For large fi the corrections can be quite large, which is reasonable since 
directly enters the Higgs boson and squark coupling vertices. They are generally larger than the 
SUSY-QCD corrections due to gluino exchanges |jlO[|- Fig. 8 shows that 6a /a is also sensitive to 
rriAO and is more sensitive for smaller tan (3. 

In Figs. 9-11, we present Sa/a for y/s = 500GeV, m^^ = 150Gel^, m^^ = 450GeV^. For these 
mass values of the stops ^/s = 500GeV is close to the peak for ti pair production and also to that 
for bi pair production at large tan/3 as shown in Fig. 6. Fig. 9 shows 6a /a as a function of n 
for tan/3 = 1.5 and Fig. 10 shows 6a /a as a function of fi for tan/3 = 0.6. We can see that for 
tan /3 < 1 the cross section for stop production is greatly suppressed. We find the cusps in the two 
figures are a threshold effect mainly coming from Fig. 5(i) when ^ 250GeV. Fig. 11 shows 
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the correction as a function of m^o for /i = AOOGeV and several values of tan From this figure 
we also see that the corrections are negative and the cross section is suppressed for tan [3 < 1. 

We then discuss a scenario with large SUSY parameters at ^/s = 2000(jel^. We take rrij^ = 
AOOGeV and m^-^ = SOOGeV in the following discussions. Figs. 12 and 13 show the ratio of the 
corrections to the tree level result for a(e+e^ —>■ titi, ^2^2) as a function of fi and m^o respectively 
when 1712 = lOOOGeV^. The cusp in Fig. 13 mainly comes from Fig. 5(c) when mH+ ^ lOOOGeV^. 
The corrections are large for small tan/3. The corrections for cr(e~^e~ ^2^2) show a singularity 
which stems from the wave function renormalization for ^2 at "^a" = 4:00GeV, where m^^ = 
+ m^o. Such singularity was also mentioned in [|1^,^]. The corrections can even reach up to 
20% in this case. 

In Fig. 14 we give the corrections to the sbottom production rates for small and large tan/3 
scenarios. When tan (3 = 30 the corrections for bi and 62 are both positive and those for 62 can be 
larger than 20%. For tan/5 = 2 the corrections tend to increase the 62 production and decrease 
the bi production. 

Finally, we will compare the contributions from three classes of diagrams, i.e., (1) the vertex 
and squark wave function corrections from Higgs bosons, (2) those from charginos and neutralinos, 
and (3) corrections to the gauge boson propagator. Especially we will show the importance of 
the third class. We find this contribution is sensitive to the mass difference between the two 
scalar top quark mass eigenstates. In Fig. 15 we give three classes of contributions as functions 
of parameter /x for = 1500^61/, tan/5 = 3, m^o = 400GeV, m2 = SOOGeV, rrif^ = SOOGeV 
and mj- = 500G'eV, SOOGeV. In fig. 16 we give the same quantities as functions of parameter 
rrif^ by fixing m^- = SOOGeV, fi = —SOOGeV and take other parameters the same as those in Fig. 
15. We find that the contribution from corrections to gauge boson propagator can be as large as 
about 7% of the total cross section. It can be seen from these figures that the contribution from 
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corrections to the gauge boson propagator is larger than those from charginos and neutrahnos in 
a large region of the parameter space and it is opposite in sign and not much smaller than the 
contributions from Higgs bosons. 



V. SUMMARY AND CONCLUSION 



In summary, we have calculated the large Yukawa coupling corrections to the diagonal stop and 
sbottom pair production in e^e~ annihilation. We include also terms of the self-energy corrections 
to gauge bosons enhanced by large masses. We discuss the corrections as functions of different 
SUSY parameters. They are found to be quite significant and are larger than the SUSY-QCD 
corrections by gluino exchanges in a large region of the MSSM parameter space. They can 



even be comparable to the conventional QCD corrections which is about 20% [|T0|- The corrections 
can be both positive and negative. We find the corrections are quite sensitive to the parameters /x, 
niAo and tan/5. They are not sensitive to the gaugino mass m2. In conclusion, when one consider 
the third generation squark production in the MSSM such corrections should not be ignored if 
precision prediction is needed. 



After we finished the work we became aware of the work of H. Eberl et al ||2^ in which a 
large part of this work had been done. However, they did not include the contributions coming 
from corrections to the gauge boson propagator. As discussed at the end of the last section this 
contribution is sizable in a large region of parameter space and should be taken into account for 
consistency. The corrections to the gauge boson propagator has been calculated in a different 
renormalization scheme in connection to the process of chargino pair production [^. The effects 
of charginos, neutralinos and Higgs bosons in the loop are not considered in that paper. Apart from 
this difference we agree with the analytical formulas given in their paper. Our numerical results 
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contain studies of effects of varying different supersymmetric parameters which are not given in 
ref ^5|]. We have checked that by taking the same parameter values as in and neglecting the 
corrections to the gauge boson propagator we obtain numerical results close to theirs. 
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Appendix A 



In this appendix we list the relevant pieces of the SUSY Lagrangian in terms of the mass 
eigenstates. We follow the conventions of ref l^^, where the full Lagrangian and the complete set 
of Feynman rules for the MSSM are given. Some abbreviations of the vertex couplings are defined 
here, which will appear in the analytic expressions in next appendix. 



= (J-^) 



^^-r-(u — At cot 3) 
2Mw 

2M, 



w 



-(/i cot 13 + A 



t , 



k % J 



Mw tan^ OwB^^ ( 6'^ + ^—^ZfzY ] + 



mf 



M\Y sin P 



mt 



2Mw sw 

= —ig2^ijk , 



^ i-Mw sin 2/3 + ^ tan /5 + ^ cot [3] Z^ Zf + ^^"^^ 



(A.l) 
(A.2) 



(A.3) 



zl'zf 

t 



ig2Dij , 



w 



Z^'Z}' 

b I 



(A.4) 



mf 



ml 



M, 



w 



16 



C-. 



= {j- 



= m 



M, 



2Mw 
2Mw 



(/X tan + 74b) 



+ 



COS/:* " " " " ■ ■ 
ciy 3 Mvi/smp J 



+ 



Mvi^ sin (3 
Zfz-A Pn 



yli rylj 



_ «5'2 



y ^ Mwcosp ^ J 



+ 



nil 



ZfZ+A Pr 



(A.5) 
(A.6) 
(A.7) 



(A.8) 



(A.9) 



(A.IO) 



(A.11) 



(A.12) 



In the above expressions, Zr, Z^ and Z'^ and Z are the mixing matrices for the two neutral 
CP-even Higgs bosons, neutraUnos and charginos respectively. 
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COS a — sm a 



sm a cos a 



1 



B 



R 



cos(a + /3), k = 1, 
— sm(a + P), k = 2, 



(A.13) 

(A. 14) 
(A.15) 



and 



tan 2a = tan 2/5 



m 



io + M| 



m 



AO 



Ml 



(A.16) 



Appendix B 



In this appendix we give some analytic results in our calculations. The vertex corrections in 
Eq. and (|3^ ) are given by 



A^ 



9l 



( m 



+ 



(47r)2 [3 
/ nit ^"^ 
\2M^ 



V2M, 



w 



(/i — At cot [3Y (C'o + 2Ci)[m|, s, m|, m^o, m^^, m 



(/i cot /3 + A)' (Co + 2Ci) [m| , m| , M|, m?^ , m?J 



+ (riafcri„fc) (Co + 2 Ci) [m|, s, ml, m%o, 171-^,7111 
+ ( Aj)^ (Co + 2 Ci) [m| , s, ml , m|, , , m^+] 
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V2M, 
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+(Ai)'(0.5 - s^)(Co + 2 Ci)[m|, s, m|, mf., m|+, ml,+] 
+ (Z^:,.)'(0.5 - sl){Co + 2 Ci)[m|, m? , m|,, M^, M^] 
-(Aj)^Cij(Co + 2 Ci)[m|, s, ml,m%+, ml mU 



{D[^fG,,{Co + 2 Ci)[m|, m|, M^, m|, m|] 



rrit 



mw sin (5 



(m^om^o — m? — 



Ci - m^Co) [m~,,s, mf., , mjo, mjo] - |5o[s, mjo, m^o 



2 _2 



{RijSij)m^omt{l - ^s^)(Co + 2Ci)[m|, s, m|, m^g, m^^ mf 



+ 4-^(4 + ^ 



3 

,2 I _2 



(m| + m'^q)Ci + m^oCo ] s, m|, m^o, , m^] + ^-Bo[s, m^, m^' 



+ (-^i/ - -Sw (Ro^ + -^i/)") mlCi[m}.,s, m\,mlo,ml, m]] 



--cos(26'h/) 

+ (^ + ^3) 



2{UijVij)m +mb{Co + 2Ci) [mf. , s, mf. , m^, m^+, m^+] 
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1 
2 



+ rrif. + m^+)Ci + m^Coj [mf.,s, mf.,m^,m^+, m^+] + -Bo[s, m^+,m^+] 
4 

{UijVij)m.^+mb{l - -s^)(Co + 2Ci)[m|, s, ml,mK, ml, ml] 



((("^1 + + ^l+'^oj K-., s, ml,ml+,ml, ml] + ^Bo[s, ml, 

+ (Vi/ - (f^i/ + "^i/)) rnlCi[ml,s, ml,m^^+,ml, ml] 



mu 



The analytic expressions for stop self energies are 



9l 



(47r)^ 



(1 - 5^^ ) 



(/i - At cot/3)^i?o[p^,r«|,,m^o] 



+(1 - b'') 



2Mi 



2Mw. 

2 

(/X cot /3 + AtfBolp^, ml , m|] 



+ {Tiak^jak)Bo[p'^,ml^,mHo] + {DjaDai)Bo\p^,ml^,mH+] + {Dj^D'^i)Bo\p^,rr 



3 + 2s? 



m? cot /? m? tan^ 



-2 sin cos ^,-5^^ ( ^^-^^ cos(2/5) - ^^^^^ + ^^^^^ ^oK±] 



12c2^ 



2m^ 



2m^ 



3 - 8s?, 



+ sing^-cosgt ^^ (5^^' (cos(2/3)^[mio] +cos(2Q;)(^[m^o] - Ao[m^o])) 



+ (RjkRik + SjkSik) (^Ao[m^o] + m^Bolp^, ml, m^o 
+ {RjkRik + SjkSik) p^Bi\p^, ml, mjo] 



+ {UjkUik + VjkVik) ( Ao[m'^-] + mlBo[p\ mi. 



m 



+ {UjkUik + VjkVik) p^Bi\p^, ml, m 



The analytic expressions for the gauge boson mass corrections are 



5Ml 5M^ N^glml 



92 



Ml 



Mi 



w 



647r2M^ 167r2Af,i 



|sin^(Q; - (5)Bqq[MI, m\o,m]jo] 



+cos^(q; - (5)Bqq[MI, m%,mlo] 
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sin^(a - /?)Boo[M^,m|-+,m^o] + cos^(q; - /?)-Boo[M^, m^+, m^o] 



+4 



-2 



+4 



Boo[Ml,ml ml 



4s^ , 3-8s2,,^i,,2' 



+ 



12 



6 b "^ly" 



Soo[M|,m? ,m 



•sw"^ _^ 3 - 4g^ ^^i^^2 



9 



12 



Ao[m 



6iJ 



'6 ""t 

,2 



{Z^^Z^^ - Zf}Z^^) -AQ[mlo] - (m^om^o + mlo)Bo[Ml mlo,m\] 



+25oo [M|, 



-m^_5o[M^,m2_,m^o]) 



-2m^om^- (Z^'Z+'^Zf Z-'')So[M^, 



(B.4) 



Vijk, Dij, D'^j, Rij, Sij, Uij and Vij in the above expressions are the vertex coupUngs defined in 
Appendix A. In the concrete calculations we only keep the higgsino sector in the R^j, S^j, Uij and 
Vij. The other two coupling constants 



G 



-Zl'Zl^ - -slyS'^ 

2 b b ^ w 



(B.5) 
(B.6) 



are the couplings of Z boson to top squark and bottom squark respectively. The relevant scalar 
functions are defined as follows 
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Ao{m') = (iTT^)" (27r/x)^-" j d^'q^q' - m')" , (B.7) 
BM.mlml) = (iny\2n^i)'-'' J d^q[{q' - ml)({q + p,)' - m?)]-\ (B.8) 

4) 



I 2 2 2 2 2\ 



(in')-\2n^i)'-'' I d^q[(q'-ml)((q + p,r-ml)((q+p2r-ml)r' , (B.9) 



in which pij = {pi — pj) 



2 



The definitions of the tensor integrals and the relevant decompositions are given below 



T,,...,Mr- ■,PN-i,mo,-- ■,mN-i) = ^^^^^ / A ^ ' ' ' , (B.IO) 

ITT J UqUi ■ ■ ■ i/jv— 1 



with the denominator factors Dq — q^ — ml, Di — {q + pi)^ — mf (i=l,- • -jN-l) and T — B,C, D ■ ■■ 
corresponding toA'" = 2,3,4---. 

B,^Pi,B,, (B.ll) 

B^,u = 9iiuBQQ + pi^i^^ii, (B.12) 

=Pi^Ci +p2^C2 = Y.Pi^^Ci, (B.13) 

i=l 

C^,u = g^ci^Coo + Pi^i^Cn + P2iJ)2uC22 + {PliJ>2^+P2^1u)Cl2 
2 

= QixvCm + ^ Pi^Pj^Cij . (B.14) 

The analytic expressions of Ao{m'^), Bo{p'^,ml,mf) and Bi{p'^,ml,ml) can be easily obtained 
and the corresponding divergences are: 

AK) (^^^-7E + ln47r) +•••, (B.15) 

So(p^ m^, m?) = - 7b + In 47r + • • •, (B.16) 

m^, ml) = ^ _ + In 47r) + • • •. (B.17) 
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Bqq can be expressed by ^40, Bi and Bq as follows 

Soo(p^ ml, ml) = ^ [Ao{ml) + 2mlBo{p^, ml, ml) + (p^ + ml - ml)Bi{p^, ml, ml) 



6 

2 

2 2 P 

+mo + mi - y 



(B.18) 



and we can extract the divergent part 

■1 



B,o{p\mi,mi) 



m^ + mo) - —p 



A-D 



-7E + ln47r + 



(B.19) 



The contribution of the top-quark loop to the self-energy of the Z gauge boson is 



a 
An 



3^8 

3 



W) 



k'^At + (/c^ + 2ml)F{k'^, rrit, mt) - 



'8s2 r2 



(At + F{k\mt, 



mt 



;b.2o) 



where 



Fik'^, wu, wu) — — I In 

^0 



x^k"^ — xk"^ + m? — ie 



ml 



(B.21) 



For k^ < ml F{k'^,m1,ml) can be expanded as a convergent power series in ^ with as 
the first term. Therefore, mlF{k'^ ,ml,ml) is not large in this region. For k"^ > mj the term 
^* ^''dZ^i^^^^ is not enhanced either. Therefore ^* ^''dZ^i ^^^ is not large for all values of k^. 
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FIGURE CAPTIONS 



FIG. 1 The tree-level Feynman diagram for the process e+e titi. 

FIG. 2 Self-energies and counter terms for internal gauge bosons. Note that each graph represents 
four different combinations. 

FIG. 3 Feynman diagrams for gauge boson mass corrections. 

FIG. 4 Feynman diagrams for self energies of squarks and their mixing at one loop order. 

FIG. 5 Corrections to vertex ■j{Z)iiii due to Higgs boson, neutralino and chargino exchanges. 

FIG. 6 The cross section a{e~^e~ — >• iiii,bibi) as a function of the collision energy ^/s for m^^ = 
150GeV, rrii^ = 450^6^, = m^o = = mCeV and (a) tan p = 1.5, (b) tan p = 30. 

FIG. 7 Corrections Sa/a as a function of /i for e+e" iiii at ■\/s = 206GeV for rrif^ = 92GeV, 
777.^2 = SBOGeV, tan/3 = 1.5 and several values of m^o and m2. 

FIG. 8 Corrections Sa/a as a function of m^o for e'^e~ iiii at ^/s — 206GeV for m^^ = 926*6^, 
rui^ = 350^6^, // = GOOGeV, = GOOGeT/ and different tan /3. 

FIG. 9 Corrections 8a ja as a function of \x for e'^e" — > t-^\ at = 500G'ey for m^^ = 1506*6^, 

= 450GeV, tan/3 = 1.5, = GOOGeV and several values of rriAo. 

FIG. 10 Corrections 5a /a as a function of n for e+e~ — > iiii at = 500G'ey for m^^ = 1506*6^, 
= 4:50GeV, tan/3 = 0.6, m2 = 600Gel/ and several values of m^o. 

FIG. 11 Corrections 5a /a as a function of uiao for e"'"e~ — > iiii at = 500GeV" for rrif^ — 
150GeV, mi^ = 450^6^, ^l = 400^6^, = eOOGeF and several values of tan (3. 
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FIG. 12 Corrections bo jo as a function of \x for e+e t\t\^t2i2 at = 20006*61^ for m^^ = 
m^GeV, mi^ = SOOGel/, tan/? = 2, m2 = lOOOGel/ and m^o = SOOGel/, QOOGel/, ISOOGel/. 

FIG. 13 Corrections bo jo as a function of m^o for e+e" ^ ^iti,t2^2 at = 2000Gel^ for 
m^^ = 400^6^, mi^ = SOOGey, // = 12006*6^, = lOOOGey and tan/3 = 1.5, 2, 30. 

FIG. 14 Corrections 5(t/(T as a function of // for e+e~ — > 6i6i,62&2 at = 2000G'ey for 
mj^ = 4006*6^, vii,^ = SOOGeV, = lOOOGey , m^o = SOOGeF, QOOGey, 1300(76^ and 
(a) tan/5 = 2, (b) tan/3 = 30. 

FIG. 15 Corrections bo jo as a function of ^ for e+e^ t\t\ due to the three classes of con- 
tributions, that is coming from, (1) corrections to the vertex and squark hues by Higgs 
bosons (2) by charginos and neutrahnos and (3) corrections to the gauge boson propagators 
respectively. The parameters arc taken to be y^ = 1500GeV", tan/3 = 3, m^o = 400GeV^, 
m2 = SOOGel/, m^- = 'imGeV and m^- = SOOGel/, SOOGel/. 

FIG. 16 Corrections bo jo as a function of m^^ for e^e~ —>■ iiii due to the three classes of 
contributions, that is coming from, (1) corrections to the vertex and squark lines by Higgs 
bosons (2) by charginos and neutralinos and (3) corrections to the gauge boson propagators 
respectively. The parameters are taken to be y^ = ISOOGeV, tan/3 — 3, m^o = 400Gel^, 
m2 = SOOGeV, // = -300^6^ and m^- = SOOGeV. 
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